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We consider the impact of colored scalars that can couple directly to matter fields on the recently 
measured h — > 77 excess. Among all possible candidates only scalar states transforming as (8, 2, 1/2) 
and (6, 3, —1/3) under the Standard Model gauge group can individually accommodate the excess 
and remain in agreement with all available data. Current experimental constraints require such 
colored states to have an order one coupling to the Standard Model Higgs and a mass below 300 GeV. 
We use the best fit values to predict the correlated effect in h — > Z'y and di-Higgs production. We 
furthermore discuss where and how these states appear in extensions of the Standard Model with 
primary focus on scenarios of matter unification. We revisit two simple 5*17(5) setups to show that 
these two full-fledged models not only accommodate a light color octet state but correlate its mass 
with observable partial proton decay lifetimes. 



I. INTRODUCTION 

Colored scalars are frequent harbingers of new physics. For example, they are inherent to any theory of matter 
unification [1]. Phenomenologically, colored scalars that can couple directly to matter fields are of particular interest. 
This feature makes them very appealing candidates for collider physics and precision flavor studies. For example, 
some have been suggested as possible explanations of the enhanced forward-backward asymmetry in tt production as 
measured at the Tevatron [2] or the (g— 2) M anomaly [3]. More recently, it has been pointed out that some could also 
contribute significantly to the unexpectedly large CP asymmetry in the decays D° — > K + K~, ir + ir~ [4]. Moreover, in 
a framework of simple one-particle extensions of the Standard Model (SM) , these states can establish a unique link 
between collider and Planck scale physics via proton decay [5]. They might even help induce tiny neutrino masses 
through loop effects [6]. 

It has been demonstrated that light colored scalars help improve unification of gauge couplings in the Grand Unified 
Theory (GUT) setting [7-16]. It could thus be possible to establish a firm connection between the colored scalar 
mass scale and observable proton decay signatures in the minimal matter unification scenarios. Moreover, the colored 
states often reside in the same representation of the GUT group as the fields responsible for the electroweak symmetry 
breaking (EWSB). This could allow one to infer the strength of Yukawa couplings of colored states to matter. These, 
on the other hand, are subject to nontrivial constraints originating from low-energy flavor phenomenology and the 
requirement of viable fcrmion masses. What one can end up with is a highly predictive class of simple models that 
connect matter stability to low-energy phenomenology [17, 18]. 

A generic feature of massive scalars in extensions of the SM is that they couple to the SM Higgs boson. Above the 
EWSB scale, one can write marginal — the so-called "Higgs portal" [19] — operators of the form <f>^$iJ^_ff, where H 
is the SM Higgs doublet, <f> is the new scalar weak (and color) multiplet and all possible color and weak contractions 
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are assumed. After EWSB this will induce corrections to the masses of $ components (</>) and couplings of the form 
(fy<j>h, where h is the physical Higgs boson. At the one-loop level, the presence of such interactions can affect Higgs 
production and decays as measured at the LHC [20]. In fact preliminary data on the Higgs- like resonance with a 
mass of nih — 125 GeV show an excess of events in the loop-induced h — > 77 decay channel with respect to the SM 
prediction at the 2tr level, while other, tree-level dominated Higgs decays to vector bosons show good consistency 
with the SM expectations [21-23]. If it were to persist at larger significance, such a deviation would point towards 
the existence of new degrees of freedom at the EW scale that couple to the Higgs boson. 

Motivated by these observations, in the present study, we focus on the contributions of colored scalars to loop- 
induced production and decay channels of the Higgs boson through interactions of the H form, where A is 
a dimensionless coupling. In particular, we study fields with direct couplings to SM matter fields that can thus play 
an interesting role in other collider and flavor observables. We single out two fields which accommodate the observed 
h — > 77 enhancement while remaining in good agreement with the other measured decay channels. Interestingly 
enough, these fields turn out to be sextets and octets of color. This precludes them from having leptoquark-like 
couplings to the ordinary matter. 

We then proceed to show how and where these particular states appear in simple scenarios of matter unification. 
We revisit gauge coupling unification in two full-fledged models based on the SU(5) gauge group [24] to demonstrate 
the connection between the lightness of these states and observable proton decay signatures. One of the models can 
also address the (<? — 2)^ anomaly, albeit through the presence of an additional light colored scalar. 

The remainder of the paper is structured as follows. In Section II we discuss a possible explanation of the h — > 77 
excess with colored scalar fields. We analyze the impact of all the relevant experimental inputs and specify the 
numerical procedure used to generate our results. We also provide predictions for h — > and di-Higgs production 
using the best-case scenarios and associated parameters. The fields that can help accommodate the h — > 77 excess 
are shown to appear naturally in matter unification scenarios in Section III. There we also correlate the lightness of 
these states with observable proton decay signatures. Finally, we conclude in Section IV. 



II. ENHANCED DI-PHOTON SIGNAL WITH COLORED SCALARS 



We study the influence of colored scalar fields coupling to the SM Higgs doublet through A<f>^ H terms on Higgs 
production and decay signatures at the LHC. If sufficiently light, such states can significantly modify the loop- induced 
gg — > h, gg — > hh, h — > jZ and h — > 77 processes, while leaving other, tree-level dominated Higgs production and 
decay channels SM-like. In particular, the partial decay width for h — > 77 is given by [20, 25], 
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where is the Fermi constant, a the fine structure constant, g w = V Ana/ sin 9 W and 9 W the Weinberg angle. Also, 
xi = m^/ (Am 2 ) for i = W,t,4>, while relevant one-loop functions are given by 



Ai{x) = -2(2x 2 + 3x + 3{2x-l)f(x)) X - 2 , 
A 1/2 (x) = 2(x + (x-l)f(x))x- 2 , 
A (x) = - (x - f(x)) x~ 2 , 
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The first and the second term in Eq. (1) are the SM one- loop contributions from the W and the top quark, respectively. 
For rrih = 125 GeV, their numerical values are A\(xw) = —8.3 and ^1/2(2;*) = 1-4. Finally, d(r,p) are the dimensions 
of the color representations of (f>, and their electric charges, while are defined via the ^(fih interaction terms 
induced after EWSB. We fix their normalization as 
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Analogously, the parton level gg — » h cross section at partonic c.m.s. energy reads [20, 25] 

Vgg^h = <J m 2 h 5(s - ml), 
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where a s is the strong coupling constant and C(r^) is the index of the color representation of (p. The only color 
representations we consider are triplet, sextet and octet as these are the ones that can contract directly with SM 
matter fields in SU(3) space to yield a singlet. We accordingly use C(3) = 1/2, C(6) = 5/2 and C(8) = 3. 

After specifying the color and EW representations, each scenario (<&) is thus parametrized in terms of the relevant 
Higgs couplings (A^) and masses (m^) of the weak components <j> of the multiplet. Different possible weak contractions 
of the \&§WH terms will in general induce different relative contributions to A^ (and corrections to m<f,). However, 
severe experimental constraints coming mainly from the p parameter require an approximate custodial symmetry to 
be active in the EW symmetric scalar potential (c.f. [26]). It turns out that in this limit, one can uniquely determine 
all the A^ in terms of a single combination of parameters in the EW symmetric Lagrangian (A). In our numerical 
results we will therefore identify A0 simply with the coupling in Eq. (6) of the maximal weight state in any given weak 
multiplet <f>, all others being thus uniquely determined. Also, in the custodial limit, </> are almost degenerate and we 
will use rritf, to denote their common mass. 



A. Numerical Procedure 



Our goal is to confront available Higgs signal strength data with possible contributions from light colored scalars. 
Parameters used to fit the data are A^ and m^. Following [21] we define individual channel signal rates, normalized 
to their respective SM values as 

( J2 a 3^h x Bh^i \ 
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where the labels j and i denote the relevant Higgs production and detection channels. We furthermore denote the 
reported experimental values and variances of /i.; by jli and af, respectively (listed in Table I). 
A global x 2 is then defined as 
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where we neglect correlations among the various terms as they are not supplied by the experimental collaborations. 
As pointed out in [22], theoretical uncertainties are only relevant for a g i ^[ Jfh , where they amount to a relative error of 
±14%. To obtain the Oi values in Eq. (7), we add this contribution in quadrature with the experimental errors (pi) 
for each observable. For each scenario (i>), we determine the minimum of the \ 2 (Xmin)> an d define the 68% (la) 
and 95% (2a) best-fit regions as solutions to \ 2 < Xmin + A% 2 , where A% 2 are set by the appropriate cumulative 
distribution function. 



B. Data 



The relevant observables, which can be significantly affected by colored scalar contributions, are the recently mea- 
sured LHC Higgs production rates in the WW*, ZZ*, 77 and 77,7 j channels [27-35] listed in Table I. In our fit we 
combine the independent measurements for each channel using weighted average. The 77JJ channel is the only one 
with a significant contribution from vector boson fusion (VBF). The expected signal strength in this channel can be 
parametrized as [21] 

0.033 a gg ^h + °Vbf ^ft-m 

^77Ji - q qoo /jSM T SM X 75SM ' 

where a gg ^h denotes the hadronic gg — > h cross section and the ratio ay^ F /a gg ^ ~ 0.078 [36] remains almost constant 
when going from 7TeV to 8TeV c.m.s. energy at the LHC. Since VBF remains SM like in our scenarios, we can also 
identify er VBF — Cvbf- All other channels are completely dominated by gg h production alone. 
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CHANNEL 


Ai ± 5"i 


REFERENCE 


pp->ZZ*->41 


1.25 ±0.55 


ATLAS7±8 [22] 


pp^ZZ*^41 


0.85 ±0.3 


CMS7±8 [22] 


pp^WW*^41 


1.4 ±0.5 


ATLAS7±8 [35] 


pp^WW*^41 


0.7 ±0.4 


CMS7±8 [22] 


pp-^77 


1.7 ±0.55 


ATLAS7±8 [22] 


pp->77 


1.5 ±0.4 


CMS7±8 [22] 


pp^77jj 


3.2 ±0.9 


CMS7 [22] 




1.6 ±0.8 


CMS8 [22] 



Table I: Data used in the analysis. 



C. Results 

We consider contributions of colored scalars listed in Table II to Higgs production and decays. The fit parameters 
are the effective coupling A^ and the colored scalar mass m$. In the SM reference scenario, Xmm/d.o.f. = 1.84. 



SU(3) x SU(2) x 17(1) 


Xmirt 


X 2 (ZZ, WW, 77, 77 jj) 


(3, 1, 1/3) 


6.8 


(0.49, 0.25, 1.61, 4.5) 


(3,1,-2/3) 


7.2 


(0.22, 0.09, 2.20, 4.7) 


(3,1,-4/3) 


6.1 


(0.31, 0.30, 2.36, 3.1) 


(3,2,1/6) 


7.1 


(0.35, 0.17, 1.91, 4.6) 


(3,2,7/6) 


6.5 


(0.20, 0.21, 2.55, 3.5) 


(3,3,-1/3) 


7.4 


(0.04, 0.01, 2.62, 4.7) 


(6,1,-1/3) 


6.7 


(0.55, 0.28, 1.48, 4.4) 


(6,1,2/3) 


6.9 


(0.44, 0.21, 1.73, 4.5) 


(6,1,-4/3) 


7.4 


(0.07, 0.03, 2.56, 4.7) 


(6,3,-1/3) 


0.7 


(0.02, 0.04, 0.12, 0.5) 


(8,2,1/2) 


1.3 


(0.03, 0.00, 0.01, 1.2) 


SM 


7.4 


(0.04, 0.01, 2.63, 4.7) 



Table II: The list of colored scalars that couple to the SM fermions at renormalizable level and corresponding Xmin from a fit 
to Higgs production and decay measurements. Last column shows contributions to Xmin from individual observables in the fit 
whereas the last row contains the SM result. 

We notice two different behaviors in the interesting range of parameters. Obviously, to enhance h — » 77, a negative 
A^ is preferred, but this tends to lower gg — > h, affecting all the measured channels. The goal then would be to have 
as small C(r^) as possible and as big d(r^)Q^ as possible, and to have a small negative A^. It turns out, however, 
that none of the considered scalars can accommodate all the data in this way. The best candidate from the list would 
be the (3,1,-4/3) and we present its x 2 plot in the (A^m^,) plane in Fig. 1. Clearly, (3,1,-4/3) can fit gg — > h 
well, but fails to enhance h — > 77. This is typical for most of the scalars listed in Table II, as can be seen from the 
last column which shows separate contributions for all observables in the fit. For the sake of the argument, we also 
plot the scalar (3, 1,8/3) example, which fits the data perfectly but does not couple directly to the SM fermions. 

Other possibility is to allow A^C^r^) to be negative and large, so that contribution of the scalar in the loop is twice 
the contribution of the top quark but with opposite sign. In this case gg —¥ h will again have the same value as in 
the SM. This particular scenario to keep gg — > h at the SM level has been observed in Refs. [37, 38]. To accomplish 
this, we need large C(r^) since A^ is constrained by perturbativity arguments. It turns out that we have two good 
candidates listed in Table II that can accomplish the task. These are the color sextet (6,3,-1/3) and color octet 
(8,2,1/2). The corresponding x 2 plots for these states are presented in Fig. 2. Note that the optimal parameter 
space in both cases is very narrow. This practically fixes the allowed A^ for a given value of and vice versa. We 
obtain Xmm/d.o.f. = 0.35 and %^^/d.o.f. = 0.63 for (6,3,-1/3) and (8,2, 1/2) scenarios, respectively. Note that in 
this best-fit region, the sextet (octet) masses above 380 GeV (340 GeV) quickly lead to non-perturbative values of the 
coupling A > V^t. On the other hand, the remaining interesting range of masses is not yet completely excluded by 
direct searches at the LHC as recently emphasized in [4]. Furthermore, current experimental analyses assume that 
such scalars are narrow and decay 100% to pairs of jets, none of which is required or predicted by our fit to the Higgs 
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Figure 1: (left) \ 2 plot for (3,1,-4/3) state where Xmin/d.o.f. = 3.05. Black strip is a region in parameter space with the 
minimum \ 2 ■ Two other regions are the la and 2a. (right) % 2 plot for (3, 1, 8/3) state where Xmin/d-o.f. = 0.35. 
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Figure 2: (left) % 2 plot for (6,3,-1/3) state contribution where Xmin/d.o.f. = 0.35. (right) \ 2 plot for (8,2,1/2) state 
contribution where Xmin/d.o.f. = 0.63. 



data (for a recent discussion on light colored resonances escaping present experimental searches c.f. [39]). 



D. Predictions for h — > Z*y and di-Higgs production 

In the SM, the h — > Zj decay is generated at the loop level in a very similar way as h — > 77. Being extremely 
suppressed, this channel has not yet been explored by the experimental collaborations. However, once it is measured, 
it will be an extremely useful test of possible explanations of the h — > 77 anomaly since, in general, these two rates 
are related. The partial decay width for h — » Z7 including colored scalars contribution is given by [40] 



64tt 4 



1 - 



^ n / -1 \ 2(1 — (8/3) sin 2 9 W ) , vsin6 w . u N ^ ■>„ , -1 \ 

wCi{x w ,y w )-\ COS W C i/2\ x t >yt)-\ ^ i^Hmwgz^d^r^Q^m^C^x^ ,y ) 
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where yi = Am1/m 2 z and the coupling of </> to Z boson is given in units of |e|, that is gz^cji = 
2(T| — sin 2 9\y)/s\n 29w- Here TjJ represents the value of the weak isospin of <f> and in Eq. (8) we sum over 
all <f> within the given weak multiplet <£>. The relevant one-loop functions are defined as 



Co(x,y) = h(x,y), 



G 



d(x, y) = 4(3 - tan 2 9 w )I 2 (x, y) + ((1 + 2x^) tan 2 9 W - (5 + 2x -1 )) Ji(x, y), 
Ci/2(x,y) = h(x,y) - h(x,y), 

where 

'■(*■») - ^ + A (/(arl) - '<»"'» + <3? <9(I " ) " » r '» ' 

= \J x~ x — 1 arcsin ir > 1. 

The SM contributions to h — > Z7 induced by the 14 7 boson and the top quark are proportional to cos 9wC\(x~^ ,yw) — 
5.8 and 2(1 — (8/3) sin 2 9w)Ci/2(x^ 1 : y t )/cos9w = —0.3, respectively. Clearly, the SM result is dominated by the W 
boson contribution. 

We calculate fj,z-y = ^h^z-y/^f-^iz-y > f° r * ne best-fit values of coupling A^ and mass for the scalars that can 
accommodate current Higgs data. For the (8, 2, 1/2) state we predict fj,%~ — (0.90 ± 0.01) throughout the interesting 
range of scalar masses and couplings. Similarly, for the (6, 3, —1/3) state we get //z 7 = (0.70±0.02), again independent 
of the scalar mass in the interesting range of parameters. We observe that both scenarios predict a mild suppression 
of the h — > Z7 decay rate and that a 20% measurement of the relevant branching fraction could start probing the 
scalar color sextet explanation of the h — > 77 excess, something possibly in reach of the 14TeV LHC with several 
hundred fb" 1 [41]. 

Another interesting related process is the di-Higgs production via gluon fusion, i.e., gg — ► hh, which is again loop 
suppressed in the SM and thus potentially sensitive to non standard contributions. We accordingly turn our attention 
to this process. The di-Higgs production can be significantly affected in the presence of light colored fields. This is 
especially true for the regime we explore where \\<f,\ is large (A^ < 0) and colored state mass is relatively small [38]. 
We accordingly evaluate [ihh = <J gg^hh/^g^f l f l for the best-fit values of coupling A^ and mass for both scalars 
that accommodate current Higgs data in a satisfactory manner. We use the parton level loop-induced top quark 
contribution towards gg — > hh amplitude [42] as well as the loop-induced colored scalar contribution [38] to evaluate 
total cross-sections crf^ hh and o g g^ hh that are relevant for the current LHC energy reach. We integrate the parton 
level gg — > hh cross sections using the LO MSTW2008 parton distribution functions [43] with fixed factorization and 
normalization scales of 2m/,. 

At the 8TeV c.m.s. energy LHC we obtain for the (8,2,1/2) and (6,3,-1/3) states an enhancement of Hhh = 
(200 ± 60) and fihh — (140 ± 40) respectively, for the best fit regions of couplings and almost independent of the 
scalar masses in their interesting range. We have checked that these values also hold for the 14TeV LHC energy. 
We observe that both scenarios predict significant enhancement of the total cross-section for di-Higgs production via 
gluon fusion with respect to the SM predictions. The signals and discovery strategies within the regime of enhanced 
di-Higgs production at LHC have been discussed extensively in Ref. [43], where the relevant rates for various final 
states can be found for both 8TeV and 14TeV c.m.s. energy LHC. 

III. COLORED SCALARS AND MATTER UNIFICATION 

Our analysis singles out two particular colored scalars — (8,2,1/2) and (6,3,-1/3) — as potential candidates that 
can consistently address the observed enhancement in h 77 decay channel if light enough. Here we want to comment 
on matter unification scenarios that predict these scalars to be in a required mass range. Note that although both 
scalars couple directly to matter they do not mediate proton decay. 

A. Color octet 

Color octet state is an appealing source of new physics. For example, it is the only scalar beside the Higgs doublet 
that can be consistently coupled to the SM quarks within the Minimal Flavor Violation framework [44]. Octet 
production at the LHC and relevant electroweak constraint on its mass and Yukawa couplings in that particular 
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context have been extensively studied [45]. More recently, there have been numerous studies of the color octet 
influences on Higgs physics in view of LHC data and potential signals [37, 38, 46-48]. Note that in all these instances 
the presence of light octet state is simply assumed. 

There are, however, viable unifying models that predict existence of a light color octet state and correlate it with 
proton decay signatures [8, 9, 14, 49]. We discuss two particular models based on the SU(5) gauge group [24] in 
what follows. The first model [49] uses one 5-dimensional (5) and one 45-dimensional (45) scalar representation to 
accommodate charged fermion masses [50] and a set of extra fermion fields from one 24-dimensional (24p) representa- 
tion [51] to accommodate neutrino masses. The phenomenology of light octet state that resides in the 45-dimensional 
scalar representation within this particular context has been addressed in Ref. [10]. The second model uses three scalar 
representations — 5, 15 and 45 — to accommodate all fermion masses [8, 9]. These two models, beside the usual matter 
representations, also use one 24-dimensional scalar representation (24) to break SU (5) down to SU(3) x SU (2) x U(l). 

In both scenarios an upper bound on octet mass is correlated with the observable partial proton decay lifetimes. In 
view of the latest experimental data on partial proton decays modes [52, 53], in particular on p — > 7r°e + , we update 
these predictions for both models to demonstrate this correlation. 

We show in Fig. 3 a viable unification for the model with one 5-dimensional and one 45-dimensional scalar repre- 
sentation and one 24-dimensional matter representation. Numerical procedure that is used to establish unification of 
gauge couplings and implement proton decay constraints is described in detail in Ref. [12]. Here, we outline the most 
important points that lead to results presented in Fig. 3. The scalar fields are denoted as 5 = ('I'd, V&t) = (1- 2, 1/2)® 
(3, 1, -1/3), 45 = (Ai, A 2 , A 3 , A 4 , A 5 , A 6 , A 7 ) = (8, 2, 1/2)® (6, 1, -1/3)® (3, 3, -1/3)® (3, 2, -7/6)9(3, 1, -1/3)® 
(3, 1, 4/3) ffi (1, 2, 1/2) and 24 = (E 8 , S 3 , £ (3 , 2) , % 2) , S 24 ) = (8, 1, 0) ffl (1, 3, 0) ffi (3, 2, -5/6) © (3, 2, 5/6) ffi (1, 1, 0). 
The extra fermions in 24 F = (p 8 , p 3 , p {3 , 2 ) , p ( g, 2 ), P24) = (8, 1,0) © (1,3,0) ffi (3,2,-5/6) ffi (3,2,5/6) ffi (1, 1,0) are 
related through a following set of mass relations [49] 

(m„ 3 + m Pa ) , . 

m P8 = mm P3 , m P(3 2) = m P{32) = — . (9) 

Here, m is a dimensionless free parameter that describes the mass splitting between masses of ps and p 3 . We accord- 
ingly present our findings in a m vs. Mqut plane, where Mqut represents the scale of gauge coupling unification. 
Mqut is maximized through numerical procedure that varies scalar and fermion masses, in accordance with mass split- 
ting constraints of Eq. (9) , in the following ranges: 200 GeV< m^ 3 , mAi , «ia 2 > m A 4 , ™a 7 , m P3 , m P8 , m P(3 2) , m P{3 2) < 
M GUT , 10 12 GeV< m,j /T ,m A 3,TO A5 < M GUT and 10 5 GeV< 77i Sg < Af GUT [51]. 

Solid lines in Fig. 3, going from top to bottom, correspond to m-Ai = 340 GeV, wa, = 500 GeV, m Al = 5TeV and 
»iAi = 50 TeV. Horizontal dashed line is due to a constraint imposed by experimental results on proton decay through 
p — > TT°e + on unification case for rn Al = 340 GeV. Note that the difference between constraints on itia 1 = 340 GeV 
case and m Al = 50 TeV case is practically negligible as it is dominated by the difference in values of the appropriate 
unified gauge coupling at the GUT scale for these two cases. 

The most important point is that all unification scenarios below the dashed line in Fig. 3 are excluded by exper- 
imental limits on p — >• 7r°e + . The proton decay signature through p — s- n°e + channel is derived assuming that the 
Yukawa matrices for matter fields are symmetric. This assumption allows for light Ag as it prevents Ag to couple 
to a quark-quark pair [54]. That, on the other hand, renders Aq innocuous as far as proton decay constraints are 
concerned. The same assumption removes dependence on unitary redefinitions of quark and lepton fields from proton 
decay operators induced through tree-level exchange of heavy gauge bosons. Note that the mass of A 6 = (3, 1,4/3) 
needs to be below 560 GeV if it is to explain the g — 2 anomaly of muon through perturbative Yukawa couplings [18]. 
We take it to be m Ae = 350 GeV to generate Fig. 3. If A 6 mass is closer to 560 GeV the allowed GUT scale would 
be slightly raised with respect to what is shown in Fig. 3. One should also worry that such a light colored state that 
resides in the same representation of SU(5) as the octet field could spoil results for the satisfactory enhancement 
of h — > 77 signal. For example, if (3,1,4/3) is to couple to the Higgs field with the same strength as the octet 
components one would need to simultaneous fit data with both fields to address the viability of such scenario. In the 
most general case one would have a situation where the triplet and the octet have different couplings to the Higgs 
boson and different masses. 

To generate results shown in Fig. 3 we update some of input parameters with regard to what is used in Ref. [12] 
to produce partial decay width for p —> 7r°e+. We use a s (mz) = 0.1184 [55], r p _ i . Tr o e + > 1.3 x 10 34 years [52, 53] and 
a = —0.0112 GeV 3 [56]. Here a is the relevant nucleon matrix element. The predicted proton lifetime for p — > 7r°e + 
due to gauge mediation is at most a factor of 5 above the current experimental limit for the m Al = 340 GeV case 
while the proton lifetime due to scalar mediation is already at the present limit. 

The GUT scale in Fig. 3 is maximized by imposing a lower bound on proton decay mediating scalars, i.e., ^t, 
A3 and A5, to be 10 12 GeV. It has been recently shown [54] that the scalar exchange dominated proton decay in 
the models with 5- and 45-dimensional scalar representations with symmetric Yukawa couplings to matter fields 
constrains the mass of ^>t through experimental data on p — >• K + v channel to be above 1.2 x 10 13 (100 GeV/us) 




Figure 3: Viable unification for the model with 5-dimensional and 45-dimensional scalar representations and an extra set of 
fermions in 24-dimensional representation in m vs. Mgut plane at the one-loop level for central values of low-energy observables. 
Solid lines, going from top to bottom, correspond to mAi = 340 GeV, m Al — 500 GeV, mA 1 = 5TeV and niA] = 50TeV. 
Horizontal dashed line is due to the constraint imposed by experimental limit on partial proton decay lifetime through p — > Tv°e + 
on rriAi = 340 GeV case. 
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Figure 4: Viable unification for the model with 5-, 15- and 45-dimensional scalar representations in mAi vs. Mgut plane at 
the one-loop level for central values of low-energy observables. Horizontal dashed line represents a limit due to the constraint 
imposed by experimental limit on partial proton decay lifetime through p — > ir°e + . 



and 1.5 x 10 11 (100 GeV/i^) GeV in the most and least conservative case, respectively. We use T p ^, K + v > 4.0 x 
10 33 years [52], while the VEVs for 5- and 45-dimensional representations — V5 and W45 — satisfy |w 5 | 2 /2 + 12|v4 5 | 2 = v 2 , 
where v(— 246 GeV) stands for the electroweak VEV. This conclusively shows that 10 12 GeV is a reliable lower bound 
on the mass of proton decay mediating scalars for one-loop unification considerations. 

We present in Fig. 4 a viable unification for the model with 5-, 15- and 45-dimensional scalar representations 
in a m Al vs. Mgut plane at the one-loop level for central values of low-energy observables. We use 200GeV< 
m S3 ,m Al ,m A2 ,m A4 ,m A6 ,m A7 ,m p(3 2) ,m P{S 2) ,m$ a ,m< s>a < M GUT , 10 12 GeV< m^ T ,m A3 ,m A5 ,m <1 > b < Af GUT , 
10 s GeV < m Es < M GUT , where 15 = ($ a ,$ b ,$ c ) = (1,3,1) ® (3,2,1/6) 8 (6,1,-2/3). In this case the pre- 
dicted proton partial lifetime for p — > 7r°e + channel due to gauge mediation is at most a factor of 26 above the current 
experimental limit for m Al = 300 GeV. The proton lifetime due to scalar mediation, on the other hand, is at the 
present limit. 

The relevant coupling of the octet to the- would-be Higgs field h originates from the following set of SU(5) contrac- 
tions: Ai5* 5 a 45^ 7 45^ and A25* 5 ,3 45" 5 45^ ( 5 . The couplings of the neutral and charged component of the octet Ai 
to the-would-be SM Higgs, under the assumption that the SM doublet primarily originates from 5-dimensional repre- 
sentation, are A^o = 2Ai + A2 and A^ = 2Ai, respectively. Any mixing between the doublets in 5- and 45-dimensional 
representation can be easily accounted for. In any case, we need to go to the limit v$ > W45. To reproduce the setup 
used in Section II where the h — > 77 excess is accounted for via custodial symmetric color octet loops one needs to 
assume that A2 is much smaller than Ai- 
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B. Color sextet 

The color sextet (6,3,-1/3) resides in 50- and 70-dimensional representations of SU(5) [57]. It has been shown 
that it can provide for gauge coupling unification within an SU(5) framework [13]. It is, however, difficult to connect 
its lightness to proton decay or correlate its Yukawa couplings with the origin of masses for matter fields without 
additional assumptions in that particular setup. This seems, instead, to require an 50(10) embedding scheme. 
Namely, the color sextet is part of 126- and 210-dimensional representations of 50(10). Both representations are 
frequently used in model building with the former one being crucial in explaining observed fermion masses. For 
example, colored scalars, which are part of 126-dimcnsional representation of 5*0(10) and are light enough to be 
accessible at the LHC, have recently been proposed in Refs. [14-16]. In fact, there already exists a viable setup with 
an intermediate-scale sextet [14]. It might thus be feasible to have very light (6, 3, —1/3) scalar that would originate 
from an 5O(10) model. We leave it to future studies. 

IV. CONCLUSIONS 

We have considered the impact of light colored scalars that can couple directly to matter fields on the recently 
observed h — > 77 excess. We find two viable scenarios where two states — (8, 2, 1/2) and (6, 3, —1/3) — can individually 
influence the excess in a constructive way and remain in excellent agreement with all available data. The colored 
states in question should have a substantial coupling to the SM Higgs of order one and a mass of order 300 GeV (or 
below) in order to explain the data. In particular, perturbativity arguments require that the sextet and octet masses 
should be below 380 GeV and 340 GeV, respectively. The best fit values for the colored scalar masses and couplings to 
the Higgs are used to generate predictions for h —> Zj and di-Higgs production. We find moderate suppression of the 
partial decay width h — » Z7 with regard to the SM value. The di-Higgs production, on the other hand, is enhanced by 
at least a factor of a hundred with respect to the SM prediction at both the 8 TeV and 14 TeV LHC. We subsequently 
study extensions of the SM where these states naturally appear with primary focus on matter unification models 
based on SU(5) gauge group. It is shown that two simple models correlate light color octet mass with observable 
proton decay. For the color octet mass of 340 GeV the predicted partial proton decay lifetime through p — > 7r°e + 
channel is a factor of 5 (25) above the current limit for the model with extra fermions (scalars) in 24-dimensional 
(15-dimensional) representation. In the model with extra fermions it is also possible to accommodate the (g — 2) p 
anomaly, albeit through a presence of an additional light colored scalar. In conclusion, the presented setup relates 
high scale matter unification and matter stability to Higgs physics via the effects of light colored scalar states. 

Note added 

After this work was completed, a study of vacuum decay constraints on colored scalars coupling to the Higgs 
appeared in [58]. The study includes the (8,2, 1/2) representation and finds that in the low energy effective theory 
comprising of the SM and the color octet weak doublet state, and for the interesting range of color octet masses as 
extracted from our analysis, vacuum stability constraints can be satisfied. 
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